function [obj_fun_rn] = risk_neutral_eval(cycle)

global U L A3 B3 MAINT_COST REPAIR_COST SHAPE2 SCALE2 SHAPE1 SCALE1

tm=(U+L)/2;

tr=(A3+B3)/2;

CM=MAINT_COST;

CR=REPAIR_COST;
age=cycle*SHAPE2*SCALE2;
state=age;
R = CR*gamcdf(state,SHAPE1,SCALE1) + CM*(1-gamcdf(state,SHAPE1,SCALE1));

sum=0;
for i=1:state
    y1=(i-1)*SHAPE2*SCALE2;
    y2=(i)*SHAPE2*SCALE2;
    sum=sum+(i)*gampdf_approx(y1,y2,SHAPE1,SCALE1);
end
L = sum+ (tr*gamcdf(state,SHAPE1,SCALE1))+(state+tm)*(1-gamcdf(state,SHAPE1,SCALE1));

% Calculate the objective function

obj_fun_rn = R/L;